METHOD FOR CALIBRATING A VOLUMETRIC FLOW METER 
HAVING AN ARRAY OF SENSORS 



CROSS-REFERENCE TO RELATED PATENT APPLICATIONS 
This application claims benefit to U.S. provisional patent application serial no. 
60/428,3 12, filed November 22, 2002 (CC-0556); U.S. provisional patent application serial no. 
60/510,765, filed October 9, 2003 (CC-0662); and U.S. provisional patent application serial no. 
60/51 1,399, filed October 13, 2003 (CC-0665), which are all hereby incorporated by reference 
in their entirety. 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The present invention relates to flow meters; and more particularly to a method and 
apparatus for calibrating a flow meter having an array of sensors arranged in relation to a pipe 
that measures a flow rate of a fluid flowing in the pipe. 

2. Description of the Related Art 

Volumetric flow measurement plays a critical role in process optimization and control 
of most industrial processes. The current industrial flow measurement market is often 
segmented into two broad technology categories: old technology and new technology flow 
meters. Old technology flow meters include flow measurement technologies that have been in 
used for greater than 70 years and include turbine meters, orifice plates and variable area flow 
meters. The new technology flow meters have emerged over the last 30-50 years and offer 
advantages over the old technologies in performance, functionality, and reliability. The major 
types of new technology flow meters include ultrasonic meters, electromagnetic flow meters, 
vortex flow meters, and coriolis flow meters. Each type has evolved to serve various aspects of 
the diverse range of applications within the industrial flow meter landscape. For example, the 
electromagnetic flow meter has emerged as the dominate type of flow meter used in the paper 
and pulp industry. 
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In particular, flow meters having an array of sensors arranged in relation to a pipe that 
measure a flow rate of a fluid flowing in the pipe are known in the art. For example, United 
States Patent No. 6,609,069, which is incorporated by reference, discloses a method and 
apparatus for determining the flow velocity of fluid within a pipe, such as an oil pipe. 

One problem with flow meters in the art is that, despite the fact that they are calibrated 
to measure a flow rate of a given fluid(s) flowing in a given pipe(s), this calibration may riot be 
correct when the flow meter is used in the field to measure different types of fluids flowing, for 
example, in different types of pipes than those in which the flow meter was originally 
calibrated to measure fluid flow. For example, the original calibration is likely to provide an 
incorrect measurement when the flow meter is used in relation to a fluid having a different 
density or viscosity than it was originally calibrated, or when the flow meter is used in relation 
to a pipe having a different wall thickness or inner diameter than it was originally calibrated, or 
when the array of sensors are spaced differently in relation to the pipe than it was originally 
calibrated. In view of this, the known flow meters are likely to have an error in their 
measurement, for which no correction is made. Moreover, the known flow meters are not 
designed to correct the original calibration based on one or more parameters that characterize 
either the array of sensors, the pipe, the fluid flowing in the pipe, or some combination thereof. 

There is a need in the prior art for a flow meter that can be calibrated in the field based 
on receiving such parameters. 

SUMMARY OF THE INVENTION 
The present invention provides a new and unique method and apparatus for calibrating a 
flow meter having an array of sensors arranged in relation to a pipe that measures a flow rate of 
a fluid flowing in the pipe. The method features the step of calibrating the flow rate using a 
calibration correction function based on one or more parameters that characterize either the 
array of sensors, the pipe, the fluid flowing in the pipe, or some combination thereof. The 
calibration correction function depends on either a ratio t/D of the pipe wall thickness (t) and 
the pipe inner diameter (D); a Reynolds number (pUD//i) that characterizes the fluid flow in the 
pipe; a ratio Ax/D of the sensor spacing (Ax) and the pipe inner diameter (D); a ratio fAx/U me as 
of usable frequencies in relation to the sensor spacing (Ax) and the raw flow rate (U me as); or 
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some combination thereof. The apparatus takes the form of a flow meter having a calibration 
correction function module performing the aforementioned functionality. 

The Reynolds number (Re), based on pipe diameter (D), characterizes many of the 
engineering properties of the flow. The Reynolds number is a non-dimensional ratio 
representing the relative importance of inertial forces to viscous forces within a flow: 

du 

viscous d u v 

Where p is the fluid density, \k is the dynamic viscosity, U is the volumetrically averaged flow 
velocity and v (= fx/p) is the kinematic viscosity. The critical Reynolds number for pipe flows, 
above which flows are considered turbulent, is -2,300. Most flows in the paper and pulp 
industry have Reynolds number ranging from one hundred thousand to several million, well 
within the turbulent regime. In addition to demarcating a boundary between laminar and 
turbulent flow regimes, the Reynolds number is a similarity parameter for pipe flows, i.e. flows 
in geometrically similar pipes with the same Reynolds number are dynamically similar. 

The flow rate being calibrated may include a volumetric flow rate (Q) or a velocity of 
flow. In operation, the volumetric flow rate (Q) is determined based on the equation: 

Q = A*U av , 

where A is a cross sectional area of the pipe's inner diameter and U av is an average flow 
velocity. 

The average flow velocity (U av ) is determined based on the equation: 

U av = the calibration correction function * U me a S , 

where U me as is a measured flow rate. 

The velocity of flow is determined by using a K-w plot. 
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The array of sensors includes an array of pressure sensors, as well as an array of strain 
or temperature sensors. 

One advantage of the present invention is that the flow meter according to the present 
invention can be calibrated in the field based on receiving the one or more parameters that 
5 characterize either the array of sensors, the pipe, the fluid flowing in the pipe, or some 
combination thereof. 

Another advantage of the present invention is that the flow meter according to the 
present invention provides a more accurate measurement after being calibrated in the field. 

Still another advantage of the present invention is that the flow meter according to the 
10 present invention can be designed and used in relation to a class of pipes, such as Sch 10 or 40, 
and calibrated in the field based on a family of calibration curves. 

In one particular application of the present invention relates to providing a new flow ' 
measurement technology well-suited for the paper and pulp industry. The present invention 
provides robust, high-accuracy, volumetric flow rate measurement for a broad range of single 
15 and multiphase flows. The present invention can be implemented with ported pressure 

transducers or with non-intrusive sensors clamped-on to existing process piping. This first- 
principles flow measurement methodology utilizes an array of sensors to listen to the unsteady 
pressure field within standard process flow lines. Sonar array processing techniques are 
employed to track the speed at which coherent structures, inherent within the turbulent pipe 
20 flow of the process fluid, convect past the sensor array. The present invention results in a 

meter performance on paper and pulp slurries ranging from 0-5% pulp. Results from a single 
phase calibration facility are also presented demonstrating 0.5% accuracy for pipes ranging 
from 3 to 16 inches in diameter. 

25 BRIEF DESCRIPTION OF THE DRAWINGS 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram of a flow meter having an array of sensors disposed axially 
along a pipe for measuring the volumetric flow of the process flow flowing in the pipe, in 
accordance with the present invention. 
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Fig. 2 is a cross-sectional view of a pipe with fluid flowing therein having an array of 
pressure sensors disposed axially along the pipe, in accordance with the present invention. 

Fig. 3 is a flow chart illustrative of the algorithm of the signal processor, in accordance 
with the present invention. 

Fig. 4 is a k-w plot showing a convective ridge illustrative of the convective flow within 
the pipe, in accordance with the present invention. 

Fig. 5 is a graph of an attenuation profile showing a dependence of a ratio of the pipe 
wall thickness (t) in relation to the pipe inner diameter (D). 

Fig. 6 is a graph showing a dependence of the Reynolds number (Re) in relation to an 
offset (%) of the flow velocity. 

Fig. 7 is a graph of the Reynolds number (Re) in relation to an offset (%) of the flow 
velocity showing a dependence of the sensor spacing and inner pipe diameter. 

Fig. 8, including Figs. 8a, 8b, 8c, show .... 

Fig. 9 is a graph of the Reynolds number (Re) in relation to an offset (%) of the flow 
velocity in the form of a family of calibration curves. 

BEST MODE FOR CARRYING OUT THE INVENTION 
Figs. 1 and 2 show a flow meter generally shown as 10 that measures the flow rate, 
including the volumetric flow rate and flow velocity, of a single phase fluid 12 (e.g., gas and 
liquid) and/or a multi-phase mixture 12 (e.g., process flow) flowing through a pipe 14 shown in 
Fig. 2. The multi-phase mixture may be a two-phase liquid/vapor mixture, a solid/vapor 
mixture or a solid/liquid mixture, gas entrained liquid or a three-phase mixture. The scope of 
the invention is not intended to be limited to the type or kind of fluid, flow, media or mixture 
being measured by the flow meter 10, or the type or kind of the pipe 14 in which the fluid, 
flow, media or mixture is flowing. 

The flow meter 10 includes a sensing device 16 comprising an array of pressure sensors 
(or transducers) 18-21 spaced axially along the outer surface 22 of the pipe 14 (Fig. 2), having 
a process flow 12 propagating therein. The pressure sensors 18-21 measure the unsteady 
pressures produced by vortical disturbances 118 (Fig. 2) within the pipe 14, which are 
indicative of the velocity of the process flow 12. The output signals (Pi(t)- P N (t)) may be 
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conditioned in one or more ways, including amplified and/or converted in an analog-to-digital 
(A/D) converter, although the scope of the invention is not intended to be limited to any such 
conditioning. 

The output signals (Pi(t) - P N (t)) of the pressure sensors 18-21 are then provided to a 
processor 24, which processes the pressure measurement data to determine the volumetric flow 
rate and flow velocity. The measurements are derive by the processor 24 by interpreting the 
unsteady pressure field within the process piping using multiple transducers/sensors displaced 
axially over about 2 diameters in length. The processor 24, including the modules and 
functionality described herein, may be implemented using hardware, software or a combination 
thereof. One software embodiment envisioned comprises, among other things, a 
microprocessor based architecture having a microprocessor, memory, input/output devices and . 
address and data bus coupling the same. A person skilled in the art would be able to implement 
the processor 24 without undue experimentation. The flow measurements can be performed 
using ported pressure transducers or clamp-on, strain-based sensors. 

The flow meter 10 measures the volumetric flow rate by determining the velocity of 
vortical disturbances or "eddies" 118 (Fig. 2) propagating through the flow 12 using the array 
of pressure sensors 18-21. Similar to that shown in U.S. Patent Application 10/007,736 filed 
November 8, 2001 (CiDRA Corporation Docket No. CC-0122A). The flow meter 10 measures 
the velocities associated with unsteady flow fields and/or pressure disturbances created by the 
vortical disturbances or "eddies" 1 18 to determine the velocity of the flow 12. The pressure 
sensors 18-21 measure the unsteady pressures Pi(t) - P N (t) created by the vortical disturbances 
1 18 as these disturbances convect within the flow 12 through the pipe 14 in a known manner. 
Therefore, the velocity of these vortical disturbances 1 18 is related to the velocity of the flow 
12 and hence the volumetric flow rate may be determined, as will be described in greater detail 
hereinafter. 

The Pressure Sensors 18-21 
In Fig. 1, each of the pressure sensors 18-21 may include a piezoelectric film sensor to 
measure the unsteady pressures of the flow 12. The piezoelectric film sensors include a 
piezoelectric material or film to generate an electrical signal proportional to the degree that the 
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material is mechanically deformed or stressed. The piezoelectric sensing element is typically 
conformed to allow complete or nearly complete circumferential measurement of induced 
strain to provide a circumferential-averaged pressure signal. The sensors can be formed from 
PVDF films, co-polymer films, or flexible PZT sensors, similar to that described in "Piezo Film 
5 Sensors Technical Manual" provided by Measurement Specialties, Inc., which is incorporated 
herein by reference. A piezoelectric film sensor that may be used for the present invention is 
part number 1-1002405-0, LDT4-028K, manufactured by Measurement Specialties, Inc. 

Piezoelectric film ("piezofilm"), like piezoelectric material, is a dynamic material that 
develops an electrical charge proportional to a change in mechanical stress. Consequently, the 

10 piezoelectric material measures the strain induced within the pipe 14 due to unsteady pressure 
variations (e.g., vortical disturbances) within the process flow 12. Strain within the pipe is 
transduced to an output voltage or current by the attached piezoelectric sensor. The 
piezoelectrical material or film may be formed of a polymer, such as polarized fluoropolymer, 
polyvinylidene fluoride (PVDF). The piezoelectric film sensors are similar to that described in 

15 U.S. Patent Application Serial No. ... (Express Mail No. EV 286 928 210 US) (CiDRA Docket 
No. CC-0676), which is incorporated herein by reference. 

The scope of the invention is not intended to be limited to the kind or type of pressure 
sensors 18-21. 

20 . The Measurement of Volumetric Flow 

To measure volumetric flow, the flow meter 10 characterizes the velocity at which the 
coherent vortical structures 118 convect past the axial array of sensors 18-21. Coherent 
structures 1 18 are an inherent feature of turbulent boundary layers present in all turbulent 
flows. Unlike conventional vortex shedding meters, no internal geometry is required to 

25 generate these structures. 

The overwhelming majority of industrial process flows involve turbulent flow 12. 
Turbulent fluctuations within the process flow govern many of the flow properties of practical 
interest including the pressure drop, heat transfer, and mixing. For engineering applications, - 
considering only the time-averaged properties of turbulent flows is often sufficient for design 

30 • purposes. For sonar based array processing flow metering technology, understanding the time- 
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averaged velocity profile in turbulent flow 12 provides a means to interpret the relationship 
between speed at which the coherent structures 118 convect and the volumetrically averaged 
flow rate. 

Turbulent pipe flows 12 are highly complex flows. Predicting the details of any 
turbulent flow is problematic, however, much is known regarding the statistical properties of 
the flow. The maximum length scale of the eddies 1 18 is set by the diameter of the pipe 14. 
These structures 118 remain coherent for several pipe diameters downstream, eventually 
breaking down into progressively smaller eddies until the energy is dissipated by viscous 
effects. 

Experimental investigations have established that eddies generated within turbulent 
boundary layers convect at roughly 80% of maximum flow velocity. For pipe flows, this 
implies that turbulent eddies will convect at approximately the volumetrically averaged flow 
velocity within the pipe 14. The precise relationship between the convection speed of the 
turbulent eddies 1 18 and the flow rate for each class of meters can be calibrated empirically as 
described below. 

In Fig. 2, the relevant flow features of the turbulent pipe flow 12 are illustrated along 
with the axial array of sensors 18-21. As shown, the time-averaged axial velocity is a function 
of radial position, from zero at the wall to a maximum at the centerline of the pipe. The flow 
12 near the wall is characterized by steep velocity gradients and transitions to relatively 
uniform core flow near the center of the pipe 14. Vortical structures 1 18 are superimposed over 
time averaged velocity profile. These coherent structures contain temporally and spatially 
random fluctuations with magnitudes typically less than 10% percent of the mean flow velocity 
and are carried along with the mean flow. Experimental investigations have established that 
the eddies 118 generated within turbulent boundary layers remain coherent for several pipe 
diameters and convect at roughly 80% of maximum flow velocity. 

From a volumetric flow measurement perspective, the volumetrically averaged flow 
velocity is of interest. The volumetrically averaged flow velocity, defined as the total 
volumetric flow rate, Q, divided by the cross sectional area of the conduit, A, is a useful, but 
arbitrarily defined property of the flow. In fact, given the velocity profile within the pipe, little 
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flow is actually moving at this speed. The precise relationship between the convection speed of 
turbulent eddies and the flow rate is determined experimentally through calibration for each. 

The Calibration Correction Module 38 

The processor 24 includes a calibration correction function module 38 for calibrating 
the flow rate using a calibration correction function based on one or more parameters that 
characterize either the array of sensors 18-21, the pipe 14, the fluid 12 flowing in the pipe 14, 
or some combination thereof. The calibration correction function depends on either a ratio t/D 
of the pipe wall thickness (t) and the pipe inner diameter (D); a Reynolds number (pUD//x) that 
characterizes the fluid flow in the pipe; a ratio Ax/D of the sensor spacing (Ax) and the pipe 
inner diameter (D); a ratio fAx/U me as of usable frequencies in relation to the sensor spacing (Ax) 
and the raw flow rate (U meas ); or some combination thereof. 

The Reynolds number, based on pipe diameter, characterizes many of the engineering 
properties of the flow. The Reynolds number is a non-dimensional ratio representing the 
relative importance of inertial forces to viscous forces within a flow: 



du 

inertial ^ dx UD 

Re = forces = ~ = 

viscous d u v 



Where p is the fluid density, fi is the dynamic viscosity, U is the volumetrically 
averaged flow velocity and v (=fi/p) is the kinematic viscosity. Pipe flows with Reynolds 
numbers exceeding a critical value, typically -2300, are turbulent. Those with Reynolds 
numbers below this value are laminar. 

The calibration correction function module 38 may be implemented using hardware, 
software or a combination thereof. One software embodiment envisioned comprises, among 
other things, a microprocessor based architecture having a microprocessor, memory, 
input/output devices and address and data bus coupling the same. A person skilled in the art 
would be able to implement the calibration correction function module 38 without undue 
experimentation. 
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The flow rate includes a volumetric flow rate (Q) and the calibration correction function 
module 38 determines the volumetric flow rate (Q) indicated as 40' based on the equation: 

Q = A*U av , 

where A is a cross sectional area of the pipe's inner diameter and U av is an average flow 
velocity. The-average flow velocity (U av ) is determined based on the equation: 

Uav = the calibration correction function * U mea s, 

where Umeas is a measured flow rate. 

The Reynolds number p\JD/fi is defined by a ratio of the fluid density (p), the 
volumetrically averaged flow velocity (U) and the pipe inner diameter (D) in relation to the 
dynamic viscosity of the fluid (jtt). 

The calibration correction function module 38 also may receive as inputs one or 
more parameters via line 38a, consistent with that discussed below. 

Fig. 3: The Basic Steps 

Fig. 3 shows the basic steps generally indicated as 46 performed by the processor 24 in 
Fig. 1, including a step 48 for performing an FFT of pressure signals Pi(t) - P N (t); a step 50 for 
determining power of the pressure signals in a K-u plane; a step 52 for determining a 
convective ridge in the K-co plane; a step 54 for calculating velocity of flow (Vc(t)) and/or 
volumetric flow (VF); and a step 62 for calibrating the velocity of flow (Vc(t) and/or 
volumetric flow (VF). The step 48 is performed by the FFT modules 30-33; the steps 50, 52, 
54 are performed by the array processor 36 and the step 62 is performed by the calibration 
correction function module 38. 

In Fig. 1, the flow meter 10 has an array of at least three acoustic pressure sensors 18- 
21, located at three locations xi,x 2 ,x 3 axially along the pipe 14. One will appreciate that the 
sensor array may include more than three pressure sensors as depicted by pressure sensor 21 at 
location x N . The pressure generated by the vortical disturbances 1 18 may be measured through 
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pressure sensors 18-21. The pressure sensors provide pressure time- varying signals 
Pi(t) 5 P2(t),P3(t),P N (t) to a signal processor 24 to known Fast Fourier Transform (FFT) logics 
30-33, respectively. The FFT logics 30-33 calculate the Fourier transform of the time-based 
input signals Pi(t) - P N (t) and provide complex frequency domain (or frequency based) signals 
5 Pi(w) ) P2(w),P 3 (co) ) Pn(cj) indicative of the frequency content of the input signals. Instead of 
FFT's, any other technique for obtaining the frequency domain characteristics of the signals 
Pi(t) - PN(t) may be used. 

The frequency signals Pi(co) - P N (w) are fed to an array processor 36, which provides a 
flow signal 40 indicative of the volumetric flow rate of the process flow 12 and a velocity 
10 signal 42 indicative of the velocity of the process flow. 

One technique of determining the convection velocity of the vortical disturbances 118 
within the process flow 12 is by characterizing the convective ridge of the vortical disturbances 
using an array of unsteady pressure sensors or other beam forming techniques, similar to that 
shown in U.S. Patent No. 6,609,069, entitled "Method and Apparatus for Determining the Flow 
15 Velocity Within a Pipe", which is incorporated herein by reference. 

The flow metering methodology uses the convection velocity of coherent structure with 
turbulent pipe flows 12 to determine the volumetric flow rate. The convection velocity of these 
eddies 1 18 is determined by applying arraying processing techniques to determine the speed at 
which the eddies convect past an axial array of unsteady pressure measurements distributed 
20 along the pipe 14, similar to that used in the radar and sonar fields. 

The array processing algorithms determine the speed of the eddies 1 18 by 
characterizing both the temporal and spatially frequency characteristics of the flow field. For a 
series of coherent eddies convecting past a fixed array of sensors, the temporal and spatial 
frequency content of pressure fluctuations are related through the following relationship: 

25 
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Here k is the wave number, defined as k=27r/X and has units of 1 /length, o> is the temporal 
frequency in rad/sec, and U con vect is the convection velocity. Thus, the shorter the wavelength 
(larger k) is, the higher the temporal frequency. 

In array processing, the spatial/temporal frequency content of time stationary sound 
fields are often displayed using "k-co plots". K-w plots are essentially three-dimensional power 
spectra in which the power of a sound field is decomposed into bins corresponding to specific 
spatial wave numbers and temporal frequencies. On a K-co plot, the power associated with a 
pressure field convecting with the flow is distributed in regions, which satisfies the dispersion 
relationship developed above. This region is termed "the convective ridge" and the slope of 
this ridge on a K-w plot indicates the convective velocity of the pressure field. This suggests 
that the convective velocity of turbulent eddies, and hence flow rate within a pipe, can be 
determined from the output of a phased array of sensor and identifying the slope of the 
convective ridge within the K-w plane. 

Fig. 4: The K-o) plot 

Fig. 4 shows an example of a KtW plot generated from a phased array of pressure 
sensors 18-21. The power contours show a well-defined convective ridge. A parametric 
optimization method was used to determine the "best" iine representing the slope of the 
convective ridge 200. For this case, a slope of 14.2 ft/sec was determined. The intermediate 
result of the optimization procedure is displayed in the insert, showing that optimized value is a 
unique and well-defined optima. 

The K-w plot shown in Fig. 4 illustrates the fundamental principle behind sonar based 
flow measure, namely that axial arrays of sensors can be used in conjunction with sonar array 
processing techniques to determine the speed at which naturally occurring turbulent eddies 
convect within a pipe. 

The Volumetric Flow Rate 
The volumetric flow rate for a given pipe is defined as follows: 

Q = A • U av 
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Where: 



A - Cross sectional area of the pipe's inner diameter 
U av = Average velocity of the flow profile within the pipe 

As described herein before, the average velocity is defined as follows: 



Uav= / 



t pUD Ax fAx 

D' ji 9 D ' Umeas 



■U, 



meas 



Where: 



Umeas = raw flow rate as determined by the algorithm 

/{} 

= calibration correction function 



_t_ 
D 

pUD 

Ax 
D 



fAx 

Umeas 



= ratio of pipe wall thickness(t ) and pipe inner diameter(D) 
: Reynolds Number 

: dependence on the ratio of sensor spacing( Ax ) and the inner diameter 
of the pipe(D) 

dependence on the useable frequencies( f ) within the array response 



Figs. 5-7 provide graphs of empirical testing that has been conducted to determine the 
correction function dependence on each of the contributing parameters. 

Fig. 5: The Pipe Wall Thickness (t) to 

Pipe Inner Diameter (D) Dependence 
Fig. 5 shows a graph indicating an attenuation profile in relation to the dependence of 
the pipe wall thickness (t) to pipe inner diameter (D). In effect, the unsteady pressure 
fluctuations created by the converting field of eddies strain the pipe at a corresponding set of 
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frequencies. These are proportional to the wavelength of the eddies. The pipe wall acts as an 
attenuator and will decrease the signal power as a function of the wall thickness and frequency. 

The Attenuation Characteristics 
The attenuation characteristics of the pipe wall change as a function of the ratio of the 
wall thickness to the inner diameter. 

Fig. 6: The Reynolds Number Dependence 
Fig. 6 shows the measured velocity of the flow is a function of the Reynolds Number. The 
graph provides a set of empirical data demonstrating this effect: The Reynolds Number effects 
the velocity profile and therefore the distribution of the velocity across the diameter of the pipe. 

Fig. 7: The Sensor Spacing to 
Pipe Inner Diameter (D) Dependence 
Fig. 7 is a graph of empirical data showing the dependence between the sensor spacing (Ax) 
and pipe inner diameter (D). In effect, the sensor spacing (Ax) defines the wavelength, or 
alternatively the wave number, of the signals to be measured through the dispersion 
relationship. This establishes the size of the vortical signals that can be measured. The 
convective velocity of this set of vorticals has a unique relationship to the average volumetric 
flow rate. 

The Dependence of Usable Frequencies 
within the Array Response 
The sensor spacing Ax defines the wavelength, or alternatively the wave number, of the 

fAx 

signals to be measured through the dispersion relationship. The ratio defines which of 

Umeas 

the measurable wavelengths should be used at a given flow rate. For the previously described 
empirical testing, the minimum and maximum values were set to 0.3 and 0.7 respectively. 
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Figure 8a shows a dimensional K-u plot; Fig. 8b shows an array gain; and Fig. 8c shows a 
non-dimensional K-w plot. 

Using the parameter fAx/U where f is the frequency in Hz (f - ^!{2 n ) where oo is the frequency 
in rad/sec) enables the k-co plot to be non-dimensionalized. This results in the ability to use the 
same region of the array gain function regardless of the velocity. This same region is 
preferably centered about the point of maximum array response (n/Ax) of the array gain 
function. (For example, per the attached figure the maximum array response occurs at k = n/Ax. 
Keeping in mind the dispersion relationship k = co/u = 2 n ff\\ the maximum array response will 
occur at different frequencies for different velocities. That is, for an array with Ax = 0.2 feet at 
3 ft/sec velocity the maximum array response will occur at 3ti/Ax = 47.1 rad/sec where at 10 
ft/sec velocity the maximum response will occur IOtt/Ax = 157 rad/sec. If the frequency is non- 
dimensionalized in this manner, the maximum array response will always occur at fAx/u = 0.5 
regardless of the velocity. The advantage of this non-dimensionalization is that regardless of 
the velocity, the section of the k- ro plot that yields the maximum array response can be 
identified and used. 

Referring to step 62 of Fig. 3, knowing the calibration function and the relationship of 
each of the factors in the calibration function, the measured velocity of the flow Uc(t) and 
measured volumetric flow (VF) signals may be calibrated knowing the actual or empirical 
relationship of each of the relationships of the calibration function to provide a more accurate 
flow velocity Uc(t) and volumetric flow (VF). . 

Fig. 9 

Fig. 9 illustrates empirical data representative of a family of calibration curves for a 
plurality of classes of flow meter of the present invention, described herein before, used to 
calibrate each respective class of meter. Each class of meter is represented or defined by the 
pipe that the flow meter is measuring. For example, one class is represented by a 4 inch ID, 
schedule 10 that is represented by the calibration curve 100. Each calibration curve using 
empirical data collected for a particular class of flow meters is defined by the equation Co + 
Ci/Re A C2 (Correction Function), wherein each class have a unique value for each of the 
coefficients CO, CI and C2. Re is the Reynolds number of the flow propagating through the 
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pipe. The Calibration Curves are defined by Percent Error or Offset versus the Reynolds 
number of the fluid, which is related to the velocity of the fluid flow as shown in Fig. 9. For 
example, V C orrected=V me asured/(Offset+l), the Offset=C 0 +Ci/RE A C 2 , and Re=0.00774 x p x V x 
ID/fi, where 0.007742=Combined conversion factor V (fVs to m/s) and ID(in to m). 

Each calibration curve includes all the factors or terms of the calibration correction 
function described hereinbefore. All meters within the same class are calibrated using the same 
calibration curve, defined by the correction function. 

As shown in Fig. 9, the measured velocity is corrected in accordance with the 
corresponding calibration curve. For example, the measured output velocity for a meter 
mounted to a 4 inch ID, 40 schedule pipe having a flow with a Reynolds number of 5.0E+05 is 
corrected for a 4.50% offset. 

Each sensor head includes 3 numbers representative of each of the three coefficients 
CO, CI and C2 of the "Correction Function". These calibration numbers are pre-programmed, 
or entered into the processor 24 via line 38a, for example, to calibrate the flow meter 10. This 
enables any sensor head 16 of a different classes to be used with any processor 24. Therefore 
the processor 24 are independent of the pipes being measured. 

The pressure sensors 18-21 of Fig. 1 described herein may be any type of pressure 
sensor, capable of measuring the unsteady (or ac or dynamic) pressures within a pipe 14, such 
as piezoelectric, optical, capacitive, resistive (e.g., Wheatstone bridge), accelerometers (or 
geophones), velocity measuring devices, displacement measuring devices, etc. If optical 
pressure sensors are used, the sensors 18-21 may be Bragg grating based pressure sensors, such 
as that described in US Patent Application, Serial No. 08/925,598, entitled "High Sensitivity 
Fiber Optic Pressure Sensor For Use In Harsh Environments", filed Sept. 8, 1997, now U.S. 
Patent 6,016,702, and in US Patent Application, Serial No. 10/224,821, entitled "Non-Intrusive 
Fiber Optic Pressure Sensor for Measuring Unsteady Pressures within a Pipe", which are 
incorporated herein by reference. In an embodiment of the present invention that utilizes fiber 
optics as the pressure sensors 14 they may be connected individually or may be multiplexed 
along one or more optical fibers using wavelength division multiplexing (WDM), time division 
multiplexing (TDM), or any other optical multiplexing techniques. 
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In certain embodiments of the present invention, a piezo-electronic pressure transducer 
may be used as one or more of the pressure sensors 15-18 and it may measure the unsteady (or 
dynamic or ac) pressure variations inside the pipe 14 by measuring the pressure levels inside of 
the pipe. These sensors may be ported within the pipe to make direct contact with the process 
flow 12. In an embodiment of the present invention, the sensors comprise pressure sensors 
manufactured by PCB Piezotronics. In one pressure sensor there are integrated circuit 
piezoelectric voltage mode-type sensors that feature built-in microelectronic amplifiers, and 
convert the high-impedance charge into a low-impedance voltage output. Specifically, a Model 
106B manufactured by PCB Piezotronics is used which is a high sensitivity, acceleration 
compensated integrated circuit piezoelectric quartz pressure sensor suitable for measuring low 
pressure acoustic phenomena in hydraulic and pneumatic systems. 

It is also within the scope of the present invention that any strain sensing technique may 
be used to measure the variations in strain in the pipe, such as highly sensitive piezoelectric, 
electronic or electric, strain gages and piezo-resistive strain gages attached to the pipe 14. 
Other strain gages include resistive foil type gages having a race track configuration similar to 
that disclosed U.S. Patent Application Serial No. 09/344,094, filed June 25, 1999, now US 
6,354,147, which is incorporated herein by reference. The invention also contemplates strain 
gages being disposed about a predetermined portion of the circumference of pipe 14. The axial 
placement of and separation distance AX,, AX 2 between the strain sensors are determined as 
described herein above. 

It is also within the scope of the present invention that any other strain sensing 
technique may be used to measure the variations in strain in the pipe, such as highly sensitive 
piezoelectric, electronic or electric, strain gages attached to or embedded in the pipe 14. 

The scope of the invention is intended to include using other types or kinds of sensors 
including ultrasonic sensors similar to that disclosed in CC-0680 (Express Mail No. EV 286 
928 152 US) and U.S. patent application serial no. 60/439,715, (CiDRA Docket No. CC-0530) 
filed January 13, 2003. The other types or kinds of sensors may sense or measure any 
parameter that converts within a flow. For example, any inhomogeneities, including 
temperature, particles, turbulence eddies, acoustic variations or disturbances. 
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Scope of the Invention 
The dimensions and/or geometries for any of the embodiments described herein are 
merely for illustrative purposes and, as such, any other dimensions and/or geometries may be 
used if desired, depending on the application, size, performance, manufacturing requirements, 
5 or other factors, in view of the teachings herein. 

It should be understood that, unless stated otherwise herein, any of the features, 
characteristics, alternatives or modifications described regarding a particular embodiment 
herein may also be applied, used, or incorporated with any other embodiment described herein. 
Also, the drawings herein are not drawn to scale. 
10 Although the invention has been described and illustrated with respect to exemplary 

embodiments thereof, the foregoing and various other additions and omissions may be made 
therein and thereto without departing from the spirit and scope of the present invention. 
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